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Abstract 

BES data on J/i/j ^ 7(7r"'"7r~7r+7r") have been analyzed into partial waves. 
We fit with resonances having J^*^ = 2++ at 1275 MeV, 0++ at 1500 MeV, 
2++ at 1565 MeV, 0++ at 1740 MeV, 2++ at 1940 MeV and 0++ at 2104 
MeV, plus a broad 0~ component. The O"*""*" resonances decay dominantly to 
aa, while 2++ resonances in the high mass region decay mainly to /2(1270)(7 
and aa; 2++ resonances from the low mass region decay dominantly to pp. 

PACS numbers: 14.40.Cs, 12.39.Mk, 13.25.Jx, 13.40.Hq 
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Typeset using REVT^ 



MARK III and DM2 reported their data on J /if) 7(7r+7r-7r+7r-) in 1986 and 1989 
respectively |1],|2|. They found that the dominant component of the spectrum below 2 GeV 
was due to J/ip '-fpp with spin-parity 0~ in the pp system, where the two p appear in 
a relative P-wave. MARK III claimed a pseudoscalar resonance at 1.55 GeV/c^. This was 
supported by DM2, and in addition they claimed that structures at 1.80 and 2.10 GeV/c^ 
were also J^'" =0^+. 

Later, E760 found three structures in 7777 with masses and widths very similar to those 
observed in J/ip ^ ^tt^ti^ti^ti^ . However, the quantum numbers ^ are forbidden to 7777. 
Both MARK III and DM2 analyzed their data in terms of J/ip —>■ jpp —>■ 77r+7r^7r+7r^, but 
they did not consider J/ip ^ ^aa 77r+7r^7r+7r^. Here a means the full vr — vr S-wave 
amplitude parametrized up to 1800 MeV. In 1995, Bugg et al. reanalyzed MARK HI 
data and added the aa decay mode to pp. Inclusion of this decay identified two of the 
peaks as / = scalar resonances decaying exclusively via aa. Those states have masses 
in the region M = 1500 and 2100 MeV/c^. An additional scalar state was required at 
M = 1750 MeV/c^, decaying dominantly to aa, but also with significant decays via pp. 

We present here an analysis of BES data on J/ip 77r+7r^7r+7r^ into partial waves. In 
this analysis, the full Monte Carlo simulation of the BES detector is used. This improves 
on the analysis of Bugg et al, which used only a simplified Monte Carlo simulation, since 
the full Monte Carlo simulation of the Mark HI detector was then no longer available. 
The full Monte Carlo simulation allows an improved study of the main background channel 
J/ip 71^71^71^71^71^ . Results are similar to those of Ref. 0, but with the significant 
addition that a 2+ amplitude is required in the high mass region around 2 GeV. 

The Beijing Spectrometer(BES) has collected 7.8 x 10^ J/ip triggers, used here. Details of 
the detector are given in Ref. ||^. We describe briefly those detector elements playing a crucial 
role in the present measurement. Tracking is provided by a 10 superlayer main drift chamber 
(MDC). Each superlayer contains four layers of sense wires measuring both the position and 
the ionization energy loss (dE/dx) of charged particles. The momentum resolution is ap/P = 
1.7%Vl + P^, where P is the momentum of charged particles in GeV/c. The resolution of 
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the dE/dx measurement is ~ ±9%, providing good tt/K separation and proton identification 
for momenta up to 600 MeV/c. An array of 48 scintillation counters surrounding the MDC 
measures the time-of-flight (TOF) of charged particles with a resolution of 330ps for hadrons. 
Outside the TOF system is an electromagnetic calorimeter made of lead sheets and streamer 
tubes and having a z positional resolution of 4 cm. The energy resolution scales bs ge/E — 
22%I\/E, where E is the energy in GeV. Outside the shower counter is a solenoidal magnet 
producing a 0.4 Tesla magnetic field. 

Candidates for the decay ^ 77r"''7r~7r"'"7r~ are selected by requiring exactly four 
charged tracks. Every track must have a good helix fit in the polar angle range —0.8 < 
cos^ < 0.8 and a transverse momentum > 60 MeV/c. A vertex is required within an 
interaction region ±20 cm longitudinally and 2 cm radially. The number of neutral clusters 
may be up to five, but only one reconstructed 7 is required in the barrel shower counter. A 
minimum energy cut of 50 MeV is imposed on the photons. Showers that can be associated 
with charged tracks are not considered. Four-constraint kinematic fits are performed to the 
final states •yyAn, •yAn and 47r. The value of x^(747r) is required to be the smallest one. The 
probability of x^(747r) is required to be larger than 5%, in order to achieve good resolution 
and to select a good photon. 

Further cuts are as follows. To remove the main background, J/ip ^ 7i^7i~^7i~7i~^7i~ , 
we required the probability Prob{x^{'yA7r)) to be larger than Prob{x^{'yjmiss4:T^)), where 
^miss means this photon is missing. If the two photons in J/ip — > j^rniss'^'^'^~'^'^'^~ a-re 
from a 7r° decay, -P^o&(x^(77miss47r)) is required to be < 1%. Next, | Umiss | = | -^miss ~ 
Pmiss |< 0.12 GeV/c^ is required in order to reject events containing more than one photon 
or containing charged kaons; here, Emiss and Pmiss are, respectively, the missing energy 
and missing momentum of all charged particles, which arc taken to be 7r"'"7r~7r"'"7r~. The 
transverse momentum of the Att system P^^ — 4 | Pmiss I ^ sin^(^^^/2) is required to be 
< 0.005 (GeV/c)^, in order to remove the background J/t/j — > n^n'^n'n'^n'; here Om-y is the 
angle between the missing momentum and the photon direction. Background J/ i]j ujti^ti~ 
events are eliminated by the cut | M^^+t^-t^o — \ > 25 MeV, in the 7r°7r"'"7r~7r"'"7r~ hypothesis 
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with only one photon detected and the 7r° associated to the missing momentum. To remove 
the small background due to J/tp ^ 'jK^K^ events, a double cut on the two tt+tt^ invariant 
masses is performed, | M^+^- — M^o \> 25 MeV. The BES data show a small signal due 
to /i(1285); this is not important for our conclusion and not of concern in this analysis, so, 
in order to remove it, an additional cut is applied to discard events within the mass region 
between 1.24 and 1.32 GeV/c^. 

The effects of the various selection cuts on the data is simulated with a full Monte 
Carlo of the BES detector; 200,000 Monte Carlo events are successfully fitted to J/ip — >■ 
7(7r+7r^7r+7r^) with Air mass below 2.4 GeV; all background reactions are similarly fitted to 
this channel. The estimated background is 31%, purely from J/ip ti^^ti^ti^ti^ti^). We 
have included this background in the amplitude analysis; it lies very close to phase space 
of J/ip ^ 7(7r+7r^7r+7r^) and has only a small effect on the fit. It includes a small signal 
(visible in Fig. ^(c) and (d) below) for 02(1320) pvr, originating from J/ip ^ a2p. When 
this 02(1320) signal is combined with a fourth pion, it does not correlate with any peak in 
the 4iT mass spectrum. Otherwise, the 4n and Sir mass distributions in the background lie 
very close to phase space and are parametrised by phase space distributions. 

Fig. |I| shows the 2tt invariant mass distribution after acceptance cuts. The acceptance 
from the Monte Carlo simulation is included into the maximum likelihood fit. Signals are 
seen due to p (or a) and /2(1270). The cr contains two components, one very broad and the 
other with a peak near the p and somewhat wider than the p. It cannot be reliably separated 
from the p as a peak. The two are distinguished by their quite different angular dependence. 
They are quite easy to distinguish in the fit, but this is not easily displayed graphically, 
because of the two combinations. The lowest peak of Fig. ||(b) is due to conversion electrons, 
and has been cut out in the fit. The Att invariant mass distribution for J/tp 77r+7r^7r+7r~ 
is shown with error bars in Fig. |^; the amplitude analysis resolves structures from 1.50 
to 2.30 GeV/c^. The dark dashed histo gram of Fig. ^ shows the 4tt mass spectrum of 
J/ip ^ 7/2(1270)0", which rises strongly above 1750 MeV. For our final fit, we use 1988 
events below a 4-7? mass of 2.4 GeV. 
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We have carried out a partial wave analysis using amplitudes constructed from Lorentz- 
invariant combinations of the 4- vectors and the photon polarization for J / ip initial states with 
helicity ±1. Cross sections are summed over photon polarizations. The relative magnitudes 
and phases of the amplitudes are determined by a maximum likelihood fit. We use i to denote 
the orbital angular momentum between the photon from J/ ip decay and the resonance in the 
production process; we use L to denote the orbital angular momentum in the decay of the 
resonance to final states pp, aa and /2(1270)(j, with L restricted < 2. Because production 
is via an electromagnetic transition, the same phase is used for amplitudes with different £ 
but otherwise the same final state. 

We now discuss resonances and possible decay modes which have been included in the fit. 
Spin-parity assignments up to J = 4 have been tried. We have tried to include vr — 7r(1300) 
or 7r — ai(1260) decay modes for 0"^ resonance. There are two decay modes of 7r(1300) 
and ai(1260): pre and avr. From the PDG 0, 7r(1300) and ai(1260) decay mainly to pvr. 
When putting this constraint into the fit, we found the contributions of the tt — 7r(1300) 
and vr — ai(1260) decay modes are negligible. The resonances and decay modes which make 
significant improvements in log likelihood are listed in Table |. Amplitudes are constructed 
in terms of the combined spin S of the decay photon and the J/ip system. 

The dominant component in the fit is 0~. We find that this may be parametrized in two 
alternative ways. The one used for the fit reported here is a simple Breit-Wigner amplitude 
of mass 1440 MeV and width 225 MeV. It is much wider than ?7(1440) and is therefore not 
to be identified with that state. The rapidly increasing phase space for decays to pp with 
L = 1 results in a very broad 0~ signal, illustrated below in Fig. H(a). An alternative more 
complicated parametrization is given in a recent coupled channel analysis of many J/ ip decay 
modes |P; it gives results almost indistinguishable from those presented here. 

The width of /2(1270) cannot be fitted with any precision. Statistics of data in the 
/2(1270) mass region are very low because the /i(1285) has been removed by a cut in 
this mass region. Hence the mass and width of the /2(1270) are constrained to the PDG 
value 0. Some 2+ component is definitely required in the mass range 1500-1700 MeV, as is 
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illustrated below in Fig. §(c). This may be attributed to /2(1565), which sits at the pp and 
ujuj threshold and can decay into 47r. In the final fit, the mass of /2(1565) has been fixed at 
1565 MeV, which is well determined by data from other sources |p [|1^] ||Tl|. The possible 
cause of our optimum mass for /2(1565) lying at 1505 MeV may be crosstalk between the 
two neighboring resonances (/o(1500) and /2(1565)). Once the masses and widths are fixed, 
both contributions are well determined. 

Comparisons with data are shown in Fig. § summed from M(47r) = 1.0 to 2.4 GeV for 
M^+^-^+^-, M^+^-, Msvr, ^pTT, cos6'^+; here is the angle of tt^ with respect to the n^n' 
pair in their rest frame, and x is the azimuthal angle between the planes of tt+tt^ pairs in 
the rest frame of the resonance X. The contributions of the various components in this fit 
are shown in Fig. ^ and Fig. ^; crosses are data and histograms the fit. 

From Fig. ^, it can be seen that 2'*' resonances in the high mass region decay mainly 
via /2(1270)(T and aa; the 2"*" resonances in the low mass region decay mainly via pp. The 
peak of /2(1565) appears at ~ 1580 MeV because of the rapidly increasing phase space for 
its 47r decay. The signals at about 1580 MeV with L=0 and L=2 are both larger than 
the combined signals in Figure ^; the reason for this is destructive interference between L=0 
and L=2. 

Branching fractions are given in Table The branching fraction for M^.^ < 3.0 GeV/c^ 
is Br{J /ip 77r+7r^7r+7r^) = (3.8 ± 0.3 ± 1.3) x 10^^. The first error is statistical and the 
second systematic. 

In order to check whether the final full fit is the best solution for BES data, four tests 
are performed. 

The first test is to check the significance level of each resonance. The 0~ contribution 
of Fig. ^ is definitely required; without it, log likelihood gets worse by 293.1. Table PTI| 
summarizes the changes in log likelihood when each resonance of the fit is dropped and 
remaining contributions are re-optimized. It indicates that the /2(1270), /o(1500), /2(1565), 
/o(1740), /2(1950) and /o(2100) should be included. 



The second test is to check the J assignment of each resonance. Table shows the 



effect of substituting alternative assignments for each resonance. It can be seen that 
making an alternative assignment for each resonance is rejected clearly. 

The third test is to scan for an additional structure. Scanning an additional structure in 
0~, 0^, 2"*", or 4+ with 2, 4, 6 and 4 more free parameters respectively indicates no further 
distinctive structures. The BES data on J/ip ^ '-jK^K^ show two possible resonance 
in the 1.7 GeV region: a scalar with mass 1781 MeV, a tensor 1696 MeV. We have therefore 
tried adding to the fit a tensor resonance in this mass region. We have found no significant 
production of spin 2 is observed in the 77r^7r^7r^7r~ channel. 

The fourth test is to fit the data in slices of mass. The data are fitted in each of 10 mass 
bins 100 MeV wide from 1350 to 2350 MeV. The comparisons of the full fit and the slice 
fit are plotted in Fig. p. The solid line presents the full fit results, while the dot stands for 
the slice analysis results. They are in fair agreement with each other. The slice analysis 
confirms strongly the presence of a tensor component in the mass region around 2000 MeV. 

The results are summarized as follows: The 0"''+ resonances decay dominantly to acr, 
while the 2"'""'" resonances in the low mass region decay dominantly to pp, and those in the 
high mass region decay to /2(1270)(J, era. 

The /o(1500) and /o(1740) are produced significantly in J/ip 77r+7r^7r+7r^; according 
to the criteria of determining the gluonic content of resonance in radiative J/ip decays 
proposed by Close et al. |T^, they may therefore contain significant glueball components. 

The BES results are approximately consistent with the results from the earlier re-analysis 
of MARK III data [4]. But the BES data favour more 2"*" in the high mass region. A possible 
cause of this difference is that a better background function for BES data is available and a 
full Monte Carlo simulation is employed in BES analysis. An important conclusion is that 
a broad 2^ resonance /2(1950) is needed around 2 GeV. This is where the 2"*" glueball may 
be expected, and J/ip radiative decays are supposed to be one of the best places to search 



for it. The /2(1950) was first observed by the WA91 group [|T4| and subsequently confirmed 
by the WA102 group in their 47r mass spectrum |]T3[; this agrees with our observation. 
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FIGURES 




m^^^- (GeV/c'} m^^ (GeV/c^) 

FIG. 1. The 2tt mass spectrum of J/i]j 77r"'"7r~7r+7r~, taking P{7rf) > ^'(Trj^) and 




FIG. 2. The An mass spectrum of J/V' ^ ^n^n vr+vr ; the dark shaded histogram is the Air 
mass spectrum of J/V' — > 7/2(1270)cr 
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FIG. 3. The comparison between data and final full fit. 
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FIG. 5. tensor component contribution. 
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FIG. 6. Numbers of events in 100 MeV slices of M(47r) (data point) from (a) 0^, (b) 2+ and (c) 
0~ and the full fit (curves); the discrepancy in phase difference (d) i;^(0^) — 1;^(0^), (e) 4'{2'^) — (pi^") 
and (f) 0(2'*') — 0(0'*') between the data and fit(radians); the value plotted is that for each slice 
minus the value from the fit. 
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TABLES 



TABLE I. Resonances fitted to the BES data and decay modes used in the final fit. 





Measured 


Used in 


final fit 


decays 


i or 25+1^2 


JP 


Resonance 


Mass 
(MeV) 


r 

(MeV) 


Mass 
(MeV) 


r 

(MeV) 






2+ 


/2(1270) 


1280ti 




1275 


185 


PP 




0" 




1440^20 




1440 


225 


PP 


e = i 


0+ 


/o(1500) 


15051^5 


14011° 


1500 


112 


GO 


i = 


2+ 


/2(1565) 


15051^° 


135ti 


1565 


131 


PP 




0+ 


/o(1740) 


1740+^^ 


12011° 


1740 


120 


aa 


e = o 


2+ 


/2(1950) 


1940ti° 


3801^,^0° 


1940 


380 


/2(1270)(7 
aa 


^D2 ^D2 


0+ 


/o(2100) 


2090^^2 


330l{°° 


2104 


215 


aa 


£ = 



TABLE II. Branching ratios for various final states in J/ij; decays, integrated up to 
M(47r) = 2400 MeV. The first error is statistical and the second systematic. 



Process 












Branching ratios 




Br{J/^p 


^7/2(1270)) 


X 5r(/2(1270) ^ 


7^^ 


TT 


"7r"'"7r ) 


(1.8 ±0.2 ±0.6) X 10" 


-4 


Br{J/^p 


^70-) 


X 5r(0~ — tt+tt" 


'it' 






(1.2 ±0.1 ±0.4) X 10- 


-3 


Br{J/%l) 


^ 7/o(1500)) 


X 5r(/o(1500) 


7^^ 




"vr+vr^) 


(3.1 ±0.2 ±1.1) X 10" 


-4 


Br{J/i; 


^ 7/2(1565)) 


X Sr(/2(1565) ^ 


7^^ 


"7r~ 


"7r"'"7r~) 


(3.2 ± 0.2 ± 1.1) X 10" 


-4 


Br{J/^p 


^ 7/o(1740)) 


X 5r(/o(1740) ^ 


7^^ 




"7r"'"7r~) 


(3.1 ±0.2 ± 1.1) X 10- 


-4 


Br{J/ij 


^ 7/2(1950)) 


X Sr(/2(1950) ^ 


7^^ 


"tt" 


"tt^tt^) 


(5.5 ± 0.3 ± 1.9) X 10" 


-4 


Br{J/ih 


^ rfo(2ioo)) 


X Z?r(/n(2100) ^ 


^^ 




-^+^-) 


(5.1 ± 0.3 ± 1.8) X 10" 


-4 
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TABLE III. Changes AS' in log likelihood function S when various components of 
the fit are dropped and remaining contributions are re-optimized. 





Resonance 


AS 


Freedom 


Significance level 


2+ 


/2(1270) 


+18.67 


5 


5.0c7 


0+ 


/o(1500) 


+21.79 


4 


>6.0cr 


2+ 


/2(1565) 


+16.32 


4 


4.8(7 


0+ 


/o(1740) 


+17.73 


4 


5.1CT 


2+ 


/2(1950) 


+37.30 


7 


>6.0cr 


0+ 


/o(2100) 


+23.58 


4 


>6.0cr 


TABLE IV. Changes AS in log likelihood function when various components are 


assi 


gned a different and all contributions are re-optimized. 




Resonance 


0- 


0+ 


2+ 4+ 


2+ 


/2(1270) 


+15.1 


+13.2 


+16.5 


0+ 


/o(1500) 


+18.0 




+9.6 +17.2 


2+ 


/2(1565) 


+15.5 


+16.0 


+11.0 


0+ 


/o(1740) 


+16.8 




+11.6 +15.7 


2+ 


/2(1950) 


+36.4 


+28.2 


+31.6 


0+ 


/o(2100) 


+13.2 




+10.0 +18.1 
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